Distinctive subdynamic features of bipartite systems 
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There are several important bipartite systems of great interest in condensed matter physics and in 
quantum information science. In condensed matter systems, the subsystems are examined tradition- 
ally by using the Green function and mean-field-like methods based on the Heisenberg representa- 
tion. In quantum information science, the subsystems are handled by composite density matrix, its 
marginals describing the subsystems, and the Kraus representation to elucidate the subsystem prop- 
erties. In this work, a relationship is first established between the two techniques which appear to 
be distinct at first sight. This will be illustrated in detail by presenting the two methods in the case 
of the celebrated exactly soluable Jaynes - Cummings model (1963) of a two-state atom interacting 
with a one-mode quantized electromagnetic field. The dynamics of this system was treated in the 
Heisenberg representation by Ackerhalt and Rzazewski (1975). We present here the corresponding 
subdynamics using Kraus representation. A relationship between the two approaches is established 
and the relative merits of the two techniques are discussed in elucidating the distinctive features 
of the subdynamics. The striking effects of interaction and entanglement are made transparent in 
two illustrative examples: the transformations that manifest in non-interacting number and spin 
representations due to interactions in their respective subspaces. 

PACS numbers: 03.65.Ca, 03.65.Ta, 03.65.Ud 



I. INTRODUCTION 



An important class of quantum many-body (QMB) 
systems involves at least two interacting species, for ex- 
ample, in condensed matter physics, electron - phonon, 
electron - photon systems are major bipartite systems 
and qubits shared by two parties A(hce) and B(ob) 
have played central roles in quantum information science 
(QIS). With the advent of recent advances in experimen- 
tal techniques, an entirely different class of problems in- 
volving time-dependent phenomena, e^., interacting bi- 
partite system of atoms and radiation [l|, Q , have opened 
up a new class of issues. The understanding of these sys- 
tems lies in formulating the respective systems in seem- 
ingly different frameworks involving quantum theory: the 
Green function formalism based on the Heisenberg repre- 
sentation in QMB condensed matter for example Q and 
the Kraus representation in QIS based on marginal den- 
sity matrix formalism [4,]. It is only recently the ideas 
of quantum entanglement and decoherence in QIS are 
making their inroads into understanding quantum phase 
transition in condensed matter [Sj- 

The purpose of the present paper is to establish a 
relationship between the two frame works, that appear 
to be distinct at first sight, by exhibiting explicitly the 
same starting point and showing how they depart in de- 
tail. Even though the total bipartite system evolves in 
time unitarily, the time evolution of the subsystems is 
quite different [4]. By introducing sub-dynamic Heisen- 
berg operators, we find a novel dynamical description of 
the subsystem time evolution. The subsystem Heisen- 
berg operators do not retain the pristine nature of the 
non-interacting particles, except at certain intervals of 



time. We illustrate this in the Jaynes- Cummings Model 
(JCM) 6] of interacting bipartite system of a two-level 
atom and a single mode radiation field. 

Quantum-limited solid-state devices have shown that 
the model of an atom in a cavity, bathed in a laser field, 
serves as a generic theoretical model for describing these 
systems. JCM serves as an exactly solvable model of 
such two interacting disparate systems - a two-level sys- 
tem and an electromagnetic field (EMF) (or more gen- 
erally, a boson field) - under certain simplifying assump- 
tions. These are the rotating wave approximation and 
weak coupling of the two-level atom with the field. JCM 
has recently been extended to strong coupling regime 0] . 
It is worth pointing out the use of JCM in several other 
situations in the recent literature: entangled systems for 
controlling symmetric qubits in trapped ions ^], linear 
optics methods to generate symmetric qubits Q, collec- 
tive atomic spin excitations [10], and stability of atomic 
clocks ill!] . 

The Heisenberg scheme in JCM was given by Ackerhalt 
et al [l2| • We present here a relationship between this and 
the QIS approach of non-unitary sub-dynamics. While 
in the QMB framework, the effects of one species on the 
other are expressed in terms of various types of renormal- 
izations, in QIS, these features are manifested differently 
via explicit subsystem operators that act in the space of 
the subsystem, details of which are elucidated by exam- 
ining atom-photon coupled system in JCM. The dynam- 
ics of the subsystems of both the photon and the atom 
exhibit consequent features of dynamical correlation. Al- 
though the sub-dynamics of the photon and spin (2-state 
atom) are generally expected to be complex for the in- 
teracting system, the Heisenberg sub-dynamic operators 



are shown here to make exphcit the physical processes, 
which deviate from the pristine dynamics. 

This paper is organized as follows: in Sec. II, the rela- 
tionship of the Heisenberg and the Kraus representations 
of subsystems is established. In Sec. Ill, detailed analysis 
on the subdynamics of JCM is given. In the final section 
IV, summary and concluding remarks are made. 



II. RELATIONSHIP OF HEISENBERG AND 
KRAUS REPRESENTATIONS OF SUBSYSTEMS 

Consider two systems A and B with the Hamiltonian 
(time-independent in this sequel) 



H = Ha + Hb + H 



AB 



(1) 



as a generic interacting composite system. Here the first 
two terms represent the Hamiltonian of the two subsys- 
tems A and B and the third term is the interaction be- 
tween them. The systems A and B are described by their 
corresponding Schrodinger operators, a(0) and 6(0) re- 
spectively. The unitary time evolution operator is given 
by 



U{t) = 



-itH 



(2) 



The Heisenberg operators an (t) and bn {t) are then given 
by 



anit) = C/^(t)a(0)(g)/st/(t) 
bH{t) = U^{t)lA<E)b{Q)U{t). 



(3) 



Here I a, Ib are respectively the unit operators in A and 
B subspaces. The time dependent Heisenberg operators 
^ are composite operators acting on the total Hilbert 
space of systems A,B. Due to unitary nature of time evo- 
lution, the pristine algebra of operators remain intact. 
Time derivatives of the Heisenberg operators obey the 
equation ih^^^ = [H,OH{t)] = HOnit) - OH{t)H. 
Coupled equations containing complicated combinations 
of operators of both subsystems result, as the Hamilto- 
nian contains terms representing the interaction between 
the two systems, thus leading to time dependent opera- 
tors in the composite space. 



A. The Heisenberg Formulation 

The formal time evolution of the total system density 
matrix p{t) in the Schrodinger representation is expressed 
in terms of the unitary evolution operator, Q, in the 
standard way: 



p{t)^U{t)p{t = 0)U\t). 



(4) 



The Heisenberg representations ^ of subsystem oper- 
ators could be realized by considering the expectation 



values of a(0), 6(0): 

Tr^Trs [a(0) ® Ib p{t)] 
Tr^Trs[/A®K0)/5(i)] 



Tr^Trs [aH{t)p{Q)] = {anit)) 
TrATTB[bH{t)pm^{bH{t)). 

(5) 



Here ( ) denotes average over the initial density matrix 
p{t — 0). This approach to realize Heisenberg repre- 
sentation for the operators, opens up a natual way to 
investigate the corresponding time-dependent subsystem 
operators - defined exclusively on the subsystem spaces - 
as will be shown in the next subsection. 



B. The Kraus Formulation 

In this formulation, one takes the initial density matrix 
p{t = 0) of the composite system to be uncorrelated i.e., 
given by a direct product of their initial density matrices: 



p(i = 0)=pA(0)®Pi3(0). 



(6) 



Substituting ([6]) in ([5]) , we deduce time evolution of op- 
erators acting solely on the subspaces as follows: 

{anit)) - Tr^Trs [a(0) (g) /i3 p(t)] 

= Tr^Trs [U{t) a(0) ® Ib U\t) pa{Q) ® /3b(0)] 

= TrA[a{t)PAm- (7) 

Here we have defined an operator that acts in the sub- 
space of A only: 

~a{t) = TrB[U\t)a{Q)®lBU{t)lA®PBm 

= TYB[aH{t)lA®PBm- (8) 

Similarly we define the operator, b{t) in the subspace B. 
We now establish the expressions for subsystem operators 
in terms of using Kraus formulation for the corresponding 
non-unitary evolution of the subsystem density matrices. 
When the initial density matrix is of the un-correlated 
form, p{t = 0) = /9yi(0)(Xi/9_B(0), the time-evolution of the 
marginal density matrices assume the following structure: 

PA{t) = TlB[U{t)pA{Q)®PBm 

= ^ V,{t)pA{0)vHt), (9) 



with ^F;(i)v^^(i) = ^A, 

i 
PB{t) = TYA[U{t)pA{Q)®PBm 

= J2 W,{t)pA{0)W^{t) 



(10) 



with Y^ W^{t)Wj{t) = Ib. 



Using the Kraus representation ^, pil)) for the time 
evolution of the subsystem density matrices we get 

Tr^Trs [a(0) (g) Ib U{t)lA ® PB(0)t/t(i)] = Tr^[a(0) pAit)] 

= TTA[a{t)pAm. 



using which we obtain, 



aW=E^/W«(0)^^W- 



(11) 



Similarly a typical time dependent operator of the sub- 
system B has the form 



^w = E^^w^(o)w^,(t). 



(12) 



Because of the non-unitary nature of the time evolution 
in the subspaces, the pristine algebra of the Schrodinger 
operators are not preserved. The dynamical aspects of 
the interactions are thus manifested in these new forms. 
It is important to note that in contrast to the Heisen- 
berg representation which makes these Schrodinger op- 
erators operate in the total space independent of the ini- 
tial state of the system, the subdynamic operators give 
more detailed information concerning the way in which 
the partner subsystem influences it. This helps us to un- 
derstand how entanglement or correlation within the full 
many body system operates in contrast to the hitherto 
known renormalization ideas. Unlike in the composite 
Heisenberg representation, time derivative of subsystem 
operators d{t), hit) in (fTTj) or (fT^ do not offer simple 
algebraic equations. 

In the next section we illustrate the above formalism 
by examining the well known JCM of interacting one- 
mode radiation field and a two-level atom. This is an 
exactly soluble model and thus offers a clear example 
where all the features derived formally in this section 
can be illustrated. 



III. SUBDYNAMICS IN JAYNES-CUMMINGS 
MODEL 

The Hamiltonian of the JCM is given by. 



Ir. 



Wo 



H = uj {a)a + —)®lA + — lR®crz+9 {a®(y+~^o) ®a^). 

(13) 
The first term in Eq. p^ is the Hamitonian of the one 
mode radiation with frequency uj ; the second term is that 
of the two-level atom with ujq as the energy difference 
between its ground and excited states, and the last term 
is the dipole interaction between the radiation and the 
atom in the rotating wave approximation. The various 
operators here have the following properties: Radiation 
field is represented by the Schrodinger creation, a^, and 
destruction, a, operators of quantized photons, obeying 
the commutation relations: 



/b, (14) 

Ir. — 2_^ l'^)('^li a^a\n)—n\n) 



[a,a'\ = iR 

oo 



The two level atom is described by the Pauli spin opera- 
tors 



1 

1 ' ' "'^ 



1 
-1 



, I A 



-i 
i 

1 
1 



(15) 



^± = 2 K±*^y), a,|T> = |T), ^.U) = -U). 



Ackerhalt et. al. [I4I express p3|) as H ~ ujN + C, in 
terms of two commuting composite operators 

N = a^a ® Ia + Ir ® <J+<y- (16) 

C — ~-{uj — Wq) Ir ® (^z + gia ® cr+a^ cr_). 

(17) 

which are constants of motion. Heisenberg represen- 
tation of photon, atom operators are then expressible in 
closed form [T^l, as will be described later for purposes 
of comparing them with the subdynamic versions of the 
same operators. 

We first express the exact solutions of JCM in the fol- 
lowing form [l3|: 

i?|0,i) = r)(o,i)|0,i), f](0i) = i(c.-c.o) = 4?^, 



H\cj>{n,s)) - r!(n,s)|</)(n,s)), 

Q,{n,s) — cj (n + -) + (3 — 2s) A„, 



An 



+ <?'("+ 1) 



1/2 



(18) 
(19) 



n = 0, l,2,...,.s = 1,2; 

10(71,1)) = cos6'„ |n+ l,i) +sin6'„ |n, t), 
10(71,2)) = -sin6l„ |n-|- 1,1) +cos6'„ |n, t), 



tan 9n = 






(20) 



n=0 



The unitary time evolution operator U{t) — e '* ^ is 
given in terms of these solutions 

Uit) = |0,i)e-'*"(°-^)(0,i| 

00 

+ E E l'^K5))e-^*"("^^'(0(n,s))|(21) 

n=0 s=l,2 

which can be expressed in terms of photon-atom corela- 
tion factors in a physically transparent manner as, 

oc 

u{t) - E{^"'"*'"'^^<-iWi^'i)<"'ii 

n=0 

+e— *("+!) [«„(t)|„,t)(„,t I (22) 

-iwn{t) \n + 1, i) (n, T I - iwn{t) \n, t) {n + l,i \] 

(23) 



Here the photon-atom correlation factors are given by, 

Vn{t) = e~*^"*sin2 6l„ + e'^"*cos2 6l„, 

Wn{t) = sin 20„ sin A„ i. (24) 



A. Heisenberg Operators in JCM 

Ackerhalt et. at. [ij] express the two basic equations 
of motion of the photon creation and the atomic flip up 
operators respectively in the forms, 






i—+uj + 2C]aH+{t) 



-g<JH+{t) (25) 
ial,(i) (26) 



Since the operator C (see (HZ])) is a constant of mo- 
tion, exact Heisenberg operator solutions could be writ- 
ten down in terms of the linear combinations of the initial 
values of the operators. These exact solutions may be 
displayed in the following form to exhibit the composite 
nature of these solutions as follows: 

a^(i) = A{t)a\t = Q) + B{t)a+{t = Q) (27) 
cfH+{t) ^ C{t)a+{t^O)+V{t)a\t^O). (28) 

Here, the coeffecients yt(i), B{t), C{t), V^t) are operator 
functions obtained by solving Eqs. (^5]) . (^5)1 : 



At) 


= 


^u. /f+e**--f_e'*^+ 


V f+-f^ 


B{t) 


= 


/ pitf+ _ itf- \ 




C{t) 


= 


,uu. p+e'*^+-f_e'*- 


V f+-f^ 


V{i) 


= 


/ pitf+ _ itf- \ 




where f± are g: 


iven by 



(29) 



f±=C± X g-^{N-lR®lA) + 



uj — uja 



Ir®Ia. 
(30) 

I 



Note that f± are constant operators as they are specified 
by C, N of m]), dni). It may be seen that at i = 0, 
A, C are unit operators, while B, T) are zero. These have 
been used to re-examine the known physical properties of 
collapse and revivals by Narozhny et. al. |14| . by working 
out these operators in the representations employed in 

([HI), dsni). 



B. 



Kraus formulation of subdynamic operators in 
JCM 



In contrast to the Heisenberg operators of photon and 
atoms acting on the composite Hilbert space, we now give 
the Kraus formulation and display the operators in their 
respective subspaces. This is achieved by employing the 
expression for the time evolution operator, (|2ip and then 
performing the procedures outlined earlier. Explicitly, 
the reduced density matrix of the atom has the form: 



PA{t) - Y. WncP\W\ 

JV=0 

WNa = {N\U\a), 



(31) 



where we have taken the initial state of the pho- 
ton to be a pure coherent state pfl(O) — |Q;)(a|, 



W) 



V2 






n=0 



|n), a 



s'"^ with 



average 



number of photons given by M — jap and the photon 

number distribution is given by, p{n) — — — . 

Similarly, the reduced density matrix of radiation is 
obtained as, 



PR{t) = Y^ {p\)s:,s,Vss^p'r 

S,Sl,S2 = T,i 

Vss. = {S\U\S2). 



Vis, 



(32) 



From this we deduce the quasi- creation operator of the 
effective photon in the radiation subspace and the spin 
operators of the atom in its subspace: 



nUi\ - 



~a\t) 



where An{t) = 



Cnit) 
Dn{t) 



e-'"^Y [I" + l)("|V^^TTA„(t) + \n+ l){n - 1| C„(i) + \n){n\ D„(t)] (33) 

n=0 
P]^ ( Vn{t)vl+^{t) + Wn(t)Wn+l{t) \l —^ 1 + P\i ( I'n (*)<_! (0 + W„ (i)u'„-l (0 a/^^— ^ j 

ip\^ {■Wn{t)Vn-l{t) ^/n ~ ■Wn^i{t)vn{t)y/n + l) (34) 

ip]l {wn{t)v*^_i{t)Vn+ 1 - Wn-iit)v^it)^/n) . 



At i = the quasi-photon operator a'l'(i = 0) correctly 
reduces to the pristine photon creation operator 



at(i==0) = ^|n+l)(n|VnTT. 



(35) 



n=0 



The factor An represents the suppression of single pho- 
ton annihilation, while C„, Dn represent 2-photon and 
no-photon annihilations respectively. The precise nature 
of the photon-QPL will depend on the choice of the pa- 
rameters and on the initial spin state. The standard com- 



J 



N{t) - J2 (Pa)s.s, Vl, (a^a) F,,,, ^ dHt)&{t) 



mutation rule, [a, a'] = Ifi, does not hold for the photon- 
subsystem operators. These relations exhibit how the 
JCM interaction modifies the pristine field operators in 
important ways, thus leading to quasi-number represen- 
tation in the radiation subsytcm. In Eq. (|33p . the coeffi- 
cients represent the effects of the photon - atom interac- 
tion in a clear fashion, reflecting the non-unitary feature 
of the Kraus transformation on the pristine photon op- 
erator, Eg. psp . Similarly, the photon number operator 
is found to be 



Y^ [ \n){n\ n + \n){n\ {p]^ wl{t) - pl^ w^i W) - ip\i I" + l)("l Wnit)vn{t) + i pl^ \n)(n + 1| w„(i)<(i)] . 

(36) 



n=0 



In Fig. (la) we display how far the operator structures 
of the photon, Eqs. (|55)l. (|55|) deviate from their pris- 
tine non-interacting versions in the subystem, for various 
values of the system parameters. In the non-interacting 
case, {a\a{t)\a) — a, {a\N{t)\a) = jap and two sets 
of curves in Fig. (la) show the effects of interaction on 
these values as a function of gt. A similar calculation 
is now presented for the atomic operators in their Pauli 
spin form. Here the initial time representations of the 
pristine operators are given by 

a+{t = 0) = I T)a I = ^1, a,(t = 0) = I T)(T I - u)a I 

(37) 
and the corresponding subdynamic spin operators are 
given by 

a+it) = e-*(|T)a|5+(f) + U)(T|5+(f) 

+ \'i){USt{t) + \i){i\St{t)) (38) 



And, 



with 



= [^-(t)]^ 

oo 

snt) = Ep(")<(^)<-i(*) 

oc 
S2{t) = y^ p{n) Wn{t)Wn-l{t) 



n=0 



a vn+T 



Yl P("X(0^»-i(0' 



(39) 



n=0 



stit) 



+i 



J2 p{n)vl_,{t)wn{t) -^. (40) 



th 


l^zit) 


= iT)(Ti^r(t) + i 

+ ITXi 1^1(0 + 1 


i)(il^ 

iXth 


5|(t) 

%{t). (41) 


Slit) 


= 


oo 

l-2y]p{n)wl{t) 

n=0 




(42) 


sm 


= - 


oo 

1-2 > ]p{n + l)wl{t) 


(43) 


sm 


= - 


-2 


oo 

i ) ] p{n)wn{t)vl{t) —= 


— = sr{t) 



n=0 



(44) 



Two highlighting features here are: (i) the pristine spin 
representations given in Eq. (|37p get transformed in spin 
space due to photon - atom interaction, and (ii) the eigen- 
values of the new subdynamic effective spin operators 
thus have different eigenvalues and off-sets from their 
pristine values of -1-1, -1, and respectively. These are 
because of the non-unitary character of the Kraus trans- 
formation of the spin operators, the eigenvalues of dz{t) 
vs. gt. are displayed in Fig. (lb) for typical represen- 
tative model parameters of JCM. These eigenvalues are 
given by {OS)z ± {DS)z, with the off-set from zero is 
given by {OS)z = {Sf{t) + 5|(t))/2, and the dispersion 

is given by (03)^ = I (IS^^SM) + |^||2 I . i^ 

the non-interacting case, these are +1 and -1, with off- 
set equal to 0. The curves in Fig. (lb) show these and, in 
addition, display the collapse and revivals. The collapses 
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FIG. 1: (Color online) (a) Magnitude of {a\a{t)\a) and 
{a\N{t)\a} vs. gt. The set of JCM parameters used are, 
Atom: Pii(O) — 1, all others zero; Photon: A'mcan = \a\'^ ~ 
10; Interaction: dash-dotted line (red): Auj/g = 7.5, dashed 
line (blue): Acu/g = 10 and solid line (green): Auj/g = 20. 
Fig. 1 (b) displays the offsets and eigenvalues of a^ (t) vs gt 

for iVnican = 10. 



mimic constant values different from ±1, depending on 
the parameters and the revivals occur at different times 
depending on the detuning. Similar results are obtained 
when other choices of parameters were made. The off- 



set being small and the dispersion different from +1,-1 
are the results of the non- unitary evolution governing the 
sub-dynamics. 

The above subdynamic features may also be inter- 
preted in a complementary way: We have already shown 
how the pristine representations of spin and radiation 
get modified due to interaction and entanglement. We 
now exhibit this in terms of the constant of motion, 
C. Upon taking the expectation value of this operator 
over the total density matrix of the system, we obtain, 
(n) + i {p\^^ - pl^) = {N{t)) + i {a,{t)) , where N{t) and 
CTz(i) are given explicitly by Eqs. (|12p and (fH)) respec- 
tively. The back action on the photon number {N{t)) is 
such that it compensates the collapse and revivals that 
occur in {az{t)). Thus the photons and the atom act in 
tandem in such a way as to maintain the constant of the 
motion. 

IV. SUMMARY 



In this paper we have established a relationship be- 
tween the traditional Heisenberg formalism for dealing 
with the dynamics of composite many-particle systems 
and the corresponding subdynamics based on Kraus rep- 
resentation. The latter brings out the effects of interac- 
tion and entanglement more directly than in the former 
description. The best known exactly soluble model of a 
bipartite system of a two-state atom interacting with a 
one mode photon field is treated in detail and the re- 
sults are llustrated in Fig. 1. This figure displays how 
the photons and the two-state atomic states represented 
as effective spin half system are modified in novel ways 
beyond the already well-known results concerning them 
such as collapse and revivals. This is clearly brought out 
in Sec. II where the JCM is described in detail in the two 
formalisms. We believe that such features are aspects of 
non-unitary sub-dynamics evolution of any multi-partite 
system when one examines how the respective pristine 
operators of the separate subsystem operators are modi- 
fied in their respective subspaces. 
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